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Abstract 

We calculate quantum fluctuations of a free scalar field in the Schwarzschild-de Sitter space- 
time, adopting the planar metric that is pertinent to the presence of a black hole in an inflationary 
universe. In a perturbation approach, doing expansion in powers of a small black hole mass, we 
obtain time evolution of the quantum fluctuations and then derive the scalar power spectrum. 
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I. INTRODUCTION 



The cosmic microwave background that we observe today is almost homogenous and 
isotropic. The background temperature in our sky is about 2.7K with a tiny fluctuation at 
a level of about 10 _5 K. This is consistent with measurements of the matter content of the 
Universe, altogether prevailing a spatially flat universe. Inflation scenario can explain the 
homogeneity, the isotropy, and the flatness of the present Universe. Moreover, quantum fluc- 
tuations of the inflaton field during inflation can give rise to primordial density fluctuations 
with a nearly scale-invariant power spectrum which is consistent with the recent WMAP 
data on cosmic microwave background anisotropies Therefore, a general assumption 
usually made in most cosmological models is that the background metric is homogenous 
and isotropic. For example, a de Sitter metric is used in the inflationary era and a flat 
Friedmann- Robertson- Walker metric is used in the subsequent hot big bang. An interesting 
notion is the recent discovery of a dominant component in the matter content, dubbed the 
dark energy, which exerts a negative pressure to drive an accelerated expansion of the Uni- 
verse. If dark energy is a form of vacuum energy, our Universe will coast to the de Sitter 
space-time or the inflating phase again in the future. 

Apparently our present Universe is not so homogeneous and isotropic because we observe 
local non-linear structures such as stars, galaxies, clusters of galaxies, and very massive 
black holes. An appropriate space-time, for example, for a massive black hole sitting in the 
accelerating Universe, would be described by the Schwarzschild metric in the vicinity of the 
black hole and by the de Sitter metric at places far from the black hole. Presently, these 
local structures are decoupled from the Hubble flow, so it suffices to use the Friedmann- 
Robertson- Walker metric to study the large scale structures of the Universe. However, in 
the early Universe the gravitational effect of a black hole to the background metric may 
be important and should be addressed. For example, the existence of a black hole or a 
distribution of black holes at the onset of inflation can invalidate the use of a homogenous 
and isotropic background metric for the calculation of de Sitter quantum fluctuations. This 
also applies to the situation when we are very near to one of these black holes that still 
exists today or not. 

The cosmic no hair conjecture infers that the inflationary universe approaches asymptot- 
ically the de Sitter space-time till the end of inflation |2j. Nevertheless, the effects of matter 
and space-time inhomogeneities to inflation should be considered. Several authors have stud- 
ied the onset of inflation under inhomogeneous initial conditions to determine whether large 
inhomogeneity during the very early Universe can prevent the Universe from entering an in- 
flationary era |3|. It was found that in some cases a large initial inhomogeneity may suppress 
the onset of inflation Q. If the inflaton field is sufficiently inhomogeneous, the wormhole 
can form before the inhomogeneity is damped by the exponential expansion jHJ . In the case 
of inhomogeneities in a dust era before inflation, some inhomogeneities can collapse into 
a black- hole space-time [sj . Furthermore, for the inhomogeneities of the space-time itself, 
energies in the form of gravitational waves can also form a black-hole space-time 0. As a 
consequence, at the onset of inflation, the distortion of the metric by these inhomogeneities 
should be taken into account. 

With these considerations in mind, in this work we will investigate the quantum fluctu- 
ations of a free massless scalar field in the Schwarzschild-de Sitter (SdS) space-time. In the 
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static coordinate system, the line element of the SdS space-time is given by 

2_ fi 2GM 2 2 \ 2 / 2GM rj2 2\ ,2 , „2 J( V2 



ds 2 = - I 1 - - #V 1 eft 2 + I 1 - - ifV 1 dr 2 + rW, (1) 

where G = M|j, M is the mass of the black hole, and is the Hubble parameter for 
inflation. Here we use the convention with c = h = 1. As is well-known, the SdS metric has 
a black hole horizon and a cosmological horizon. In the static coordinates ([T]) an observer 
can only receive a signal inside or just right on the cosmological horizon. This static metric 
is insufficient for our purpose because in the cosmological setting we aim at studying the 
temporal evolution of a Fourier mode of the scalar quantum fluctuations that crosses the 
cosmological horizon during inflation. Therefore, we will instead use the planar coordinates 
for the SdS metric j8|, which is given by 

ds 2 = -f{r, r)dT 2 + h{r, r){dr 2 + rW), (2) 

where dr = a _1 (r)<it is the conformal time and a(t) = e . In Eq. ([2]), the / and h functions 
are given by 



f(r,r)=a 2 (r) 
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with the cosmic scale factor a(r) = —1/(Ht). Note that at large times the planar coordinates 
behave like a de Sitter expansion. Also, we can cast it to a co-moving frame that is conformal 
to the Minkowskian space-time with a black hole in it. 

It is well-known that a black hole evaporates into Hawking radiation which leads to a 
mass loss of the black hole and to its eventual disappearance |9j. For a Schwarzschild black 
hole with mass M, the evaporation time is given by 

t ev = 5120ttG 2 M 3 . (4) 

Therefore, our present consideration requires the condition that the evaporation time scale 
of the black hole is longer than the time scale of inflation, i.e., Ht ev > 1. This gives the 
lower bound on the mass of the black hole for a given inflation scale: 

> 3.96 x 1(T 2 I —j- ) . (5) 
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In the next section, we will review the scalar quantum fluctuations in the de Sitter metric. 
Then, we will introduce a perturbation method to expand the planar metric in powers of 
the black hole mass M to find an approximate solution of the scalar equation. Section IIHI 
contains the numerical results of the first-order solutions. Section IIVI is our conclusion. 



II. PERTURBATION APPROACH TO THE SCHWARZSCHILD-DE SITTER 
SPACE-TIME 

A. Classical solution 

Consider a massless scalar field <fi which satisfies the Klein-Gordon equation in the SdS 
space-time, 

9, (v=£<rM = o. (6) 
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Since the space is spherically symmetric, one can expand as 

^x) = Y,Vi(r,r)Y lm (OA)- (7) 
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If one further writes <fi(r, r) in a spectral form, 



oo 

2, 



(fi(r,r)= / dkk ji(kr)(p k i(r), (8) 
Jo 

then we will have 

p oo 

4>(x) = / dky^ipkimjx), (pkim(x) = k 2 ji(kr)ip H (T)Yi m (6,(j)). (9) 

J ° lm 



We define the spectral function of the fluctuations of the field (f>(x) as 

S kl (r,r) = (21 + l)jf(kr)P kl (r), P kl (r) = ^m(t)\ 2 . (10) 

Our next task is to calculate the function ip k i(r) in Eq. ( fill in the SdS planar metric fl2]). 
The Klein-Gordon equation ([6]) beomes 




where we have separated out the angular part of <f)(x) in Eq. ([7]). There is no exact solution 
to this equation. We therefore adopt a perturbative approach assuming that the quantity, 

e = GMH, (12) 

is a small parameter. But the condition (jHJ) implies that 

2 
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However, this shows that the smallness of e can be easily satisfied for any reasonable inflation 
scale. Let us rewrite Eq. (TTTj) as 

<9 r y?z - -<9 T y?; - o r (pi - -o r (pi H — 2 — (pi = I 1 - - 1 d T ifi — — II — I <9 T ^ 



7t^ t ( t2 yy) + 7k 5r (14) 



We then expand the functions / and h in powers of e as defined in Eq. ( IT2|) and write 

(°) i (!) i (2) , 

<Pi = <Pi +<Pi +<Pi +■■■ ( 15 ) 

as an expansion in orders of e. It is straightforward to show that the right-hand side of 
Eq. (Ill]) becomes 
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1. Zeroth order 



From Eqs. f|T4|) . f lT5|) . and f lT6|) . the zeroth order corresponds to the de Sitter case with 



To solve this equation, we take the Bessel transform, 



(17) 



^ (0) (r,r)= / dkeUkr^'ir) . 
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Then we have 



and the solution is 
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If we take the boundary conditions: 



then we will have 



and 
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and C 2 = 0, 
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(21) 



(22) 



(23) 



As r -> 0, P^\t) — > H 2 /(4:7T 2 ). This result matches the well-known scale-invariant power 
spectrum of de Sitter quantum fluctuations 13 • 



2. First order 



With the perturbative expansion for ipi in Eq. ( JT5]) . the Klein-Gordon equation in Eq. ((T 
can be solved perturbatively. The first order (pp (r, r) then satisfies 



(1) V„M , 
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where the source term is given by 
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To solve this inhomogeneous equation, we use the Green's function G(r, r; r', r') which sat- 
isfies the equation, 



(26) 



Using the completeness property of the spherical Bessel functions, 



dkk z 
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Eq. f[2"6"j) becomes 



POO 

Gi(r,T;r',r') = / dkk 2 g k (T,T')ji(kr)ji(kr') 
Jo 



2 2 
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T 7T 



(27) 



(28) 



(29) 



For the retarded Green's function, g k = for r' > r > 7*, where r» denotes an initial time 
when the source begins to operate. For > r > r', 



9k(r, t') 



2r' 2 k 3 
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2 2 

-(-kr')^H ( P(-kT f )(~kT^Hf\~kr) . 
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With this retarded Green's function, the first order ip^' (r, r) can be expressed as 



(30) 
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Hence, we find that 
(1) 2zeiJ 



dk'k r2 (k'i) I dr'r'ji{kr')ji{k'r') / dr' e - ik ' T ' x 
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(32) 

It is useful to rewrite <p$( T ) as 



(33) 



The integral over r' can be performed, 



dr r ji(kr')ji(k r 
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and after some rescalings the coefficients ocu and fiki can be simplified to 



-2iT{l + 1) 

w+ww 



dk'k' l+ ^F[ I + + --h 



1 



rl 



'2' ' 2' 



kr 



dT , e -ik r e ir ( r , + ^ 



kr, 



+ 



dk'k 



'u'-i+h 



3 f 



At, 



(35) 



and 
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The power spectrum in this order is then given by 



kri 



l + eA«(r 
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where we have introduced the dimensionless quantity, 
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As r -> 0, from Eq. $22 



Ag ) (0) = 2Re[aH(0)-/3 fc z(0)] 



From Eqs. ( 1351) and ( 1361) . we find that 
(1W _ 8T(Z + 1) r 
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(36) 



(37) 



(38) 
(39) 



(40) 



B. Quantization 

A unique mode function ipki m (x) can be obtained once an appropriate vacuum is chosen. 
By using this mode function the scalar field can be quantized in the standard manner, 



[x 



roc 

/ dk^\a Mm (p k i m (x) + al lm (p* klm (x) , (41) 

J ° lm 1 J 
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with the commutation relations: 



[0-klm, a>k'Vm'\ 
it 



S(k - k')8i V 8 n 



0. 



[aklm,a' k 'l'm'l = - K )°ll'°mm'- (42) 

We have the delta function in k because there is no coupling between different /c-modes. 
We will show this explicitly later by a perturbation approach. The delta functions in I and 
m stem from the rotational invariance about the central black hole. The vacuum state is 
defined as 

a k i m \0) = 0. (43) 

Now ipkim(%) is the mode function, so it should also satisfy the Klein-Gordon equation. 
Since the space is spherically symmetric, one can write 

l -Pkim(x) = <fki(r, r)Yi m (9, <p). (44) 

But (pki{ r , T ) cannot in general be separated as a product of functions with only one variable 
like k 2 ji(kr)(pki{ T ) a s we have done in the classical solution. However, ip k i(r,r) can still be 
obtained perturbatively. The Klein-Gordon equation for the mode function is the same as 
Eq. (ITT]) , given by 



Vfn 9r {\l^ dr ^ r,T ^ + ^7f^ dr ( r2 ^^(r,r)) -^4^ 



<p kl {r,r) = 0. (45) 



The two-point correlation function is then given by 
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where 7 is the separation angle between the two points. As x' — > x, we have 
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In terms of <f k i(r, T ), the spectral function of the fluctuations of the quantum field <p(x) can 
be defined as 

S kl (r,r) = (2l + l)jf(kr)P kl (r,r), P kl (r,r) = - \ n . \ip kl (r,r)\ 2 . (48) 



Perturbatively, 
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1. Zeroth order 



Following the same steps in Eqs. (JHJ), (115]) . and (1161) . to the zeroth order we have as 
before 



%<P$(r,T) - ~d T( p%>(r,T) - d 2 ^>(r,r) - -d r ^^r) + '^L^'far) = 0, (52) 
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and the general solution is found to be 



(53) 



where <p^i( T ) is given by Eq. ( |20l) . The boundary conditions ( 12T1) indeed correspond to the 
choice of the Bunch-Davies vacuum that selects the mode function: 
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and hence the zeroth order power spectrum in Eq. (|50|) is 
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2. First order 



To the next order, we have 
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The retarded Green's function necessary to solve this equation is the same as before. Then, 
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This is in the form of 
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As the perturbative formalism in quantum mechanics, the eigenf unctions in the first order 
consist of zeroth order eigenfunctions of different energies. 

Now we calculate the leading correction (l5"Tj) to the power spectrum. Substituting Eq. (I53"j) 
in Eq. (159]) . we have 
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Here we obtain from Eqs. ( )60l) and ( ISTT) that 
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(64) 



where the remaining integral over r' can be evaluated using the formula (|34|) . This quantity 
diverges as — > k'. However, we find that the integral in Eq. ( 163|) is finite and weakly 
dependent on r. In the limit of r — > 0, it can be approximated by 
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III. NUMERICAL RESULTS 



Eqs. (138]) and f )62|) are the main results of the present paper. However, they are compli- 
cated integrals and are not illuminating. Therefore, we perform numerical calculations of 
the black-hole corrections to the power spectrum. Assume that inflation begins at the initial 
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time t = and a(t = 0) = 1. Then, the initial conformal time is = —1/H and the final 
conformal time is r = —e~ Ht /H, where Ht is the e-folding number of inflation. It is useful 
to note that the Fourier mode with momentum k crosses out the cosmological horizon at 
time t = —1/k, or equivalently, when the e-folding number is Ht = \n(k/H). 

In Fig. [U we show the time evolution of the first-order contribution to the quantum 
fluctuations from Eq. (1371) . Actually, it is more convenient to plot the normalized power 
spectrum Ajy (r) in Eq. (|38[) against the e-folding number Ht. We have chosen the angular 
momentum I = 2, 22, 40 as examples. For each I, different k modes are shown. From all the 
plots, we can see that for a k-mode Aj^ (t) oscillates when the mode is still sub- horizon. 
Once the mode crosses out the horizon, A^(r) stop oscillating and gradually approaches a 
constant value. This behavior can be easily explained by Eq. f!25p . where the source term 
dies off as (p^ goes super-horizon and then gets frozen. In Fig. [21 we plot the asymptotic 
values Ajp(0) in Eq. QM) for 1 <l < 50 and 1 < k/H < 50. The figure shows that Ajp(0) is 

suppressed in low-/ and low- A; regions and the magnitude of |A^(0)| is of order one. Thus, 
the first-order contribution, compared to the zero-order de Sitter power spectrum, is roughly 
downsized by the expansion parameter e. 

For the quantum case, we plot the asymptotic values Ajy (r, 0) in Eq. (I6"3"j) for Hr = 
0.01, 0.1 against I and k/H for 1 < I < 50 and 1 < k/H < 50. As shown in Fig. [31 the 
general trend is similar to the classical case in Fig. [2) This is indeed an explicit example 
that quantum fluctuations generated during inflation behave like classical waves. We have 
also calculated A^ (r, 0) for larger values of Hr, which becomes fluctuating but in general 
the value of the amplitude is getting smaller. It is expected because the farther the black 
hole is, the lesser is its effect to the perturbation. 

IV. CONCLUSIONS 

We have presented a perturbation method to compute the effect of the presence of a black 
hole in the de Sitter space to the quantum fluctuations of a free massless scalar field. The 
method is valid as long as the expansion parameter e = GMH < 1, i.e., the product of the 
mass of the black hole and the de Sitter Hubble parameter has to be smaller than the Planck 
mass squared. The calculation can be easily applied to a vector field or a gravitational wave. 
Here the first-order contribution is computed and the results are given in the assumption 
that the black hole is located at the origin of the coordinates. Higher-order corrections can 
be worked out perturbatively though complicated. It would be interesting to consider the 
effect due to a distribution of black holes in the de Sitter space. In fact, the perturbation 
can be in a form of cosmological defects such as monopoles, cosmic strings, or domain walls. 

Let us briefly discuss some cosmological implications of the results that we have obtained 
in this work. Assume that inflation lasts for about 60 e-folds. Then, the wavelength of the 
Fourier mode with k/H = 1 is about the size of the present Universe. If the location of 
the black hole that exists during inflation is near the Earth, the suppressed power of the 
first-order correction to the de Sitter inflaton fluctuations in low I and low k regions may 
result in a blue-tilted density power spectrum on large angular scales. This in turn gives 
rise to a suppression of the large-scale cosmic microwave background anisotropy that may 
be relevant to the observed low quadrupole in the WMAP cosmic microwave background 
anisotropy data [ll[ . A detailed calculation of the effect to the cosmic microwave background 
anisotropy is underway, including the case that the black hole locates somewhere else in the 
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Universe. 

If inflation lasts for a longer time, then the wavelength of the Fourier mode that cor- 
responds to the size of the present Universe will be given by a larger value of k/H. In 
Fig. HI we plot the asymptotic values A^(0) against I for k/H = 1, 10, 100, and 250, which 
correspond to the inflation with 60, 62.3, 64.6, and 65.5 e-folds, respectively. As expected, 
the longer is the inflation duration the lesser pronounced are the effects of the black hole to 
large-scale or low-/ observations. This can also be seen in the quantum case with (r, 0) 
in Fig. [3j 

It is worth noting that the present work gives a realization of the general discussions 
in Ref. jl2j about the potentially observable effects of a small violation of translational 
invariance during inflation, as characterized by the presence of a preferred point, line, or 
plane. This violation may induce derivations from pure statistical isotropy of cosmologi- 
cal perturbations, thus leaving anomalous imprints on the cosmic microwave background 
anisotropy [l~2"j ]. 
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FIG. 1: Time evolution of the first-order normalized power spectra A^ j (r) in Eq. (pj) for I = 
2,22,40. For each case, we have selected some /c-modes. The time is expressed in terms of the 
e-folding number Ht. 
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FIG. 2: Two views of the three-dimensional plot of the asymptotic values AL (0) in Eq. (|39|) 
against I and k/H for 1 < I < 50 and 1 < k/H < 50. 
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Hr = 0.1 




FIG. 3: Three-dimensional plots of the asymptotic values AL (r, 0) in Eq. (I63p for ffr = 0.01, 0.1 
against Z and k/H for 1 < Z < 50 and 1 < k/H < 50. 
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FIG. 4: Asymptotic values A^ ; (0) against I for k/H = 1, 10, 100, and 250. 
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